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ABSTRACT 


SMITH, DEWEY M. A Comparison of Experimental Heat Transfer Coef- 
ficients in a Noaale With Analytical Predictions From Bartz’s 
Methods for Various Combustion Chamber Pressures in a r:;!»lid 
Propellant Rocket Motor. (Under the direction of JOHN NOBLE 
PERKINS.) 

The experimental heat transfer coefficients measured in the 
nozzle of a small solid propellan Tiotor are compared to the predic- 
tions from D. R. Bartz r s Nusselt i iber correlation equation and his 
technique of solving the boundary layer momentum and energy equations 
simultaneously for the heat transfer coefficient* The propellant was 
a composite of ammonium perchlorate and polybutadiene acrylic acid 
and the average motor chamber pressures were 220 , 4l0, and 742 psia. 

The nozzle was made with a steel casing and a ETA graphite throat 
insert. Measurement locations were at a local to throat area ratio 
of I .785 in the convergent section, at the throat and at a local to 
throat area ratio of 2.369 in the divergent section for the chamber 
pressures of 220 and 4l0 psia. The measurements were made at a 
local to throat area ratio of 2.067 in the convergent section, at the 
throat, and at a .local to throat area ratio of 3*764 in the divergent 
section for the chamber pressure of 742 psia. Test measurements 
consisted of temperature responses of five thermocouples mounted on a 
line perpendicular to the heated surface of the nozzle at each of 
the measurement locations. The heating rates and the corresponding 
convective heat transfer coefficients were determined by using the 
thermocouple data as input to a finite difference heat balance program. 


The results of the comparisons showed the experimental data from the 
convergent region and throat to be consistently lower than the predic- 
tions made using Barts f s two techniques. It was found that the experi- 
mental data in the divergent sections could be correlated by 
evaluating the skin friction coefficient at the free stream temperature 
and using it in the simultaneous solution. 
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INTRODUCTION 


During recent years, attempts have been made to analytically 
predict heat transfer rates in the combustion chamber and along the 
nozzle wall of solid propellant rocket motors. The prediction tech- 
niques were developed by researchers in the liquid propellant rocket 
engine field and they have been fairly successful in defining the 
heating loads for this type engine. 

The two most widely used methods for predicting the heat transfer 
rates were set forth by D. R. Bartz. The first method is a Nusselt 
number correlation equation which was developed before the advent of 
high speed computers and is still used today for rapid estimations. 

The second method solves the boundary layer momentum and energy equa- 
tions for the heat transfer coefficient and requires the use of a 
computer . 

These techniques are now being applied to solid propellant rocket 
motor nozzles. The experimental data on heat transfer rates for this 
type motor are not extensive and therefore it is difficult to determine 
the applicability of Bartz 1 s techniques to solid propellant motors. 

Most of the nozzles on solid propellant motors are not externally 
cooled and their design is based on the materials in the nozzle wall 
being able to absorb the heat transferred from the exhaust gases. 

During the times immediately after ignition, severe temperature 
gradients are set up through the nozzle wall and this condition must 
be taken into consideration in the design. Accurate prediction of the 
heating loads is desirable so that the nozzle may be designed 
efficiently. 
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To provide additional experimental data, the heating rates in a 
converging - diverging nozzle on a small solid propellant rocket motor 
were measured at a point in the convergent region, at the throat, and 
at a point in the divergent region. The average chamber pressures were 
220, 4l0, and 7^2 psia which represent the range of chamber pressures 
used in full scale motors. The heating rates were determined from 
the temperature responses of five thermocouples mounted on a line 
perpendicular to the heated surface of the nozzle at each of the 
measurement locations. The experimental heat transfer coefficients 
were then compared with predictions from the two techniques of Bartz. 
This thesis presents these comparisons and the intention is to provide 
some basis for applying the Bartz techniques to solid propellant rocket 


motor nozzles. 
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REVIEW OF THE LITERATURE 

Bartz (1+) has documented in detail the developments of both the 
analytical and the experimental research in rocket nozzle heat transfer . 
Published in 1965 , it provides a convenient catalog of the most recent 
experimental results and techniques and was referred to often in the 
performance of the work for this thesis. 

The initial prediction method of rocket nozzle heat transfer was 

based on the turbulent pipe flow heat transfer correlation equations 

of McAdams (13) * This method assumed the flow in the nozzle to be 

fully developed with each point on the nozzle contour assumed to be 

preceded by a long pipe. Sibulkin and Bartz (15, 2) were the first 

to treat the nozzle heat transfer problem with a boundary layer 

approach by making use of the Integral momentum and energy equations 

{ Sibulkin f s was an incompressible flow analysis). These treatments 

were an improvement over the McAdams equations since the flow in 

rocket nozzles is not usually fully developed. The main difference 

In Bartz’ s initial boundary layer analysis and the one used today is 

the method of solving the boundary layer equations. The original 

analysis was done before the advent of high speed computers. Bartz 

pointed out the equations were interdependent upon the ratios of wall 

temperature to stagnation temperature , =- , and the temperature to 

i o 

velocity boundary layer thickness* ^ . However, to simplify the 
mechanics of solution of the equations, he assumed initial values for 
these ratios and solved the two equations separately. The boundary 
layer equations were reduced to linear ordinary differential equations 
with variable coefficients. 
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The other basic assumptions in the analysis were: 

(1) 1/7 power law profiles of both velocity and the difference 
between stagnation temperature and wall temperature in the turbulent 
boundary layer, 

(2) the local skin-friction coefficients along the nozzle are 
the same as those on a flat plate for the same boundary-layer thick- 
ness, and 

(3) Reynolds analogy between momentum transfer and heat transfer 

applies for the nozzle boundary-layer flow. Elliot et al. (8) 

developed a computer program in 1963 that solves the boundary-layer 

equations simultaneously by an iterative method which allows the 

ratios, =■£ and — , to vary along the nozzle wall* 

To o 

A Nusselt number correlation equation was also developed by 
Bartz (3) for the purpose of making calculations of the local heat- 
transfer coefficients by hand. This method was based on the solution 
for the heat-transfer coefficient from the original boundary-layer 
analysis which showed the local coefficient to be strongly dependent 
upon the local mass flow rate. The boundary- layer solutions also 
showed that the local diameter should be used as the characteristic 
length in the correlation equation. The proportionality constant was 
obtained by matching the heat- transfer coefficients at the throat of 
a particular nozzle with the coefficient that was calculated from the 
boundary- layer analysis. This correlation equation is still used 
today for rapid calculations of the local heat -transfer coefficient. 
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Numerous experiments have been performed to determine the accuracy 
of the methods. These experiments, for the most part, have involved 
heat transfer measurements in nozzles with heated air or various liquid 
propellant exhaust gases as the working fluid- The heated air experi- 
ments are unique in that the inlet conditions to the nozzle can be 
controlled thus eliminating such combustion effects as secondary flows, 
oscillations in pressure, or free stream turbulence that can occur in 
rocket motors- Of particular interest are the data of Back, et al. (l) 
which are from heated air experiments with stagnation pressures ranging 
from 30 to 250 psia and stagnation temperatures over the 1000 -2000° R 
range. The nozzle was a 30-degree half -angle convergent, 15-degree 
half-angle divergent nozzle which is similar to the nozzles used in the 
set of tests reported on here. 

The boundary layer analysis predicted the heat transfer coefficient 
accurately throughout the nozzle whereas the correlation equation was 
approximately 50 percent high for chamber pressures in the range 75 - 
250 psia. One other interesting point was the indication that the 
boundary layer turbulence seemed to decay back toward transition of 
laminar flow near the throat in the lower stagnation pressure tests 
and was reflected in the heat transfer results. 

The data obtained by Kolozsi (10), using air at stagnation pres- 
sures of 225 and 370 psia and stagnation temperatures of 1100 - 1200° R 
in a convergent-divergent nozzle, indicated the correlation equation 
predicted coefficients too high (^5 percent high at the throat) and 
the boundary layer equations were accurate throughout the nozzle. 

For tin! and Ehlers (9) found that both methods predicted the heat 
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transfer coefficients accurately in a Rao-design divergent section 
nozzle using air at a stagnation pressure of 300 psia and a stagnation 
temperature of l600° R. Their results also showed that two-dimensional 
flow must be considered In this type nozzle instead of one- dimensional 
flow that can be used in most convergent-divergent nozzles, 

Welsh and Witte (19) used a ~ hydrazine liquid propellant 

rocket engine to gather heat transfer data and compared it with only 
the correlation equation for stagnation pressures between 80 - 290 psia. 
They found the predictions to be considerably lower in the convergent 
region when compared with the experimentally determined coefficient, 
from 80 percent above to 4 5 percent below in the throat region, with 
the best correlation in the divergent section- They theorized that 
the effects of combustion in the vicinity of the nozzle inlet 
influenced the flow in the convergent and throat region. Convcrgent- 
£ divergent nozzle configurations were used in these tests v 4 th varying 

contraction ratios • 

■Ml Witte and Harper (20) used the same liquid propellant engine used 

by Welsh and Witte (19) with nozzles over an extended range of contrac- 
tion and supersonic area ratios. They had the same general results 
^ as in (19) . They also concluded that the wide variation of data in 

the throat region of (19) was due to the -transitional tendencies of 
the turbulent boundary layer here. This tendency was felt to be 
caused by the acceleration of the How. 

Lee (12) obtained experimental heat transfer data from a solid 

. * r *. 

| propellant motor with an uncooled molybdenum nozzle that was found to 

agree with Bartz's correlation equation. The data could be correlated 
; also by assuming the skin friction coefficient to be dependent on the 
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momentum thickness and numerically integrating the boundary layer 
momentum equation, Brinsmade and Besmon ( 5 ) conducted tests with a 
solid propellant motor at stagnation pressures between l6o - 500 psia 
and a stagnation temperature of ^900° R and found that the data at the 
throat could be correlated by using laminar heat transfer equations. 
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THEORY 

The integral momentum and energy equations of the turbulent 
boundary layer can be derived by writing momentum and energy balances 
on a control volume in the vicinity of the wall where there are viscous 
effects in a real fluid. This was the method used by Bartz (4). The 
derivation is based on the definition;! of displacement, momentum, and 
energy thicknesses as the deficiencies in mass, momentum, and energy 
caused by friction and heat transfer. These deficiencies are deter- 
mined by comparing real flow with potential flow near the wall where 
the mass flow rates are made equal in the control volume for these two 
types of flow (see Fig- 1). 

The basic definitions are: 

(1) Displacement thickness 

s * ■ - % -/ o 5r (i - §)<* <» 

which is the difference in thickness in the two control vuxumes in 
order to h»ve the mass flow rates equal. 

(2) Momentum thickness 

9 • J!' r <2> 

which is the thickness of potential, flow which has a momentum flux 
that is equal to the difference between the potential and real flow 
momentum fluxee for the same ness flux. 
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(3) Energy thickness 


-- 

r ou 


t -n 

-1 o W\ - 

3- ” m m 1 ^3f 

O * W/ 


which is the thickness of potential flow that has an enthalpy flux 
equal to the difference between the enthalpy fluxes of the potential 
and real flows for the same mass flux- In forming these definitions* 
the time-mean flow density, pu, has been represented by the product of 
the mean values* p and u* and ignoring the cross- correlation terms - 
This can be done on the assumption the correlation terms cancel out 
when integrated over the boundary layer as suggested by Shapiro (l4) . 

Using these definitions in the momentum and energy balances on 
the control volumes in the potential and real flows* the following 
equations are formed: 

(l) Integral Momentum Equation 


d0 C f 9 1+ I d0 1 d(pU) 1 dr 

~ T ~ U dx pUdx rdx 


where C„ is defined as 
f 


Cf = 


(2) Integral Energy' Equation 


/ aw - w \ j 1 d(pU) . 1 dr 

" C h ( T 0 - T w J " ® T 


where C is defined as 


pUCfT 
p\ aw 


1 
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The basic assumptions used by Bartz in solving these equations 
are as follows: 

(1) The flow is axi symmetric and steady, the forces acting on 
the gas are the pressure gradient and skin friction at the wall, and 
the boundary layer is small compared to the distance from the axis of 
symmetry. 

(2) The flow through the nozzle is reversible and adiabatic with 
the change in total enthalpy of the gas due to the heat flux to the 




V*' 

■ 

... » 
n* :V.~' 


/ ■ r t f ^ : ' 

1 «| 'X-jt. . 


wall . 

(3) The gas is perfect, has a constant Prandtl number, and its 
viscosity is related to the gas temperature raised to a power. 

(*0 The skin-friction coefficient and the Stanton number are 
the same as they would be on a flat plate at the same free- stream 
conditions, wall temperature, and momentum thickness. 

(5) The Stanton number for unequal momentum and energy thicknesses 
is that for equal thicknesses multiplied by , where (n) is an 

"interaction component/ 1 The Stanton number for equal momentum and 
energy thicknesses is related to the skin- friction coefficient by 


von Harman* s form of Reynold ! s analogy 


'h 


' f ' (l-Pr + l ^ 6 
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( 8 ) 


(6) Heat transfer has either no effect on the skin- friction 
coefficient, and is the same as for adiabatic flow, or has an 

effect and the C f is the seme for adiabatic incompressible flow with 
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the density and viscosity evaluated at the arithmetic mean of the wall 
and free stream static temperatures » 

(7) The boundary layer velocity and temperature distributions 
are l/7-power profiles . 

Values for the skin friction coefficient are taken from data of 
adiabatic flow over flat plates correlated by Coles (7) . Coles found 
the data could be represented by one curve of versus C^g where 

the low speed value* C^* is related to the actual C^ by 



where T 0 is a temperature within the boundary layer which is found 



and e is the exponent in the viscosity relationship* \i ~ T 21 * The 
Hg* is related to B g by 


C fa R 0 


C f R 0 



( 11 ) 


A film temperature correction can be made by evaluating the gas 
properties p and p,* at the arithmetic mean temperature of the free 
stream temperature and the wall temperature. When this is done the 
relationship between and is 



( 12 ) 
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The term, -g- , in the momentum equation can be evaluated from the 
integrals which define 5* and 0, equations (l) and (2)* by using the 
assumed l/7-power profiles for the velocity and temperature ratio in 
the boundary layer. The limits on the integrals for these thicknesses 
as well as for the momentum thickness , 0, are expressed in terms of 5 
and A, which a^e not necessarily equal. Therefore the boundary layer 
momentum and energy equations are solved iteratively for 9 and 0 
in order to determine the proper ^ as well as the proper skin 
friction coefficient. 

The mass rate of flow per unit area through the nozzle can be 
determined from one dimensional isentropic relationships. By 
specifying the nozzle contour, wall temperature, and chamber condi- 
tions, the heat transfer coefficient at any point along the nozzle 
wall can be found from solving the boundary layer equations for the 
Stanton number, C : 




( 13 ) 


In the solutions for h developed by Bartz (2), it was found 

S 

that the heat transfer coefficient was a strong function of the mass 
flow rate per unit area, pU. From this, Bartz developed a nondimen- 
sional equation in the form 

ha * c (%.) a ("Pry m 

for the determination of h . The analysis by Bartz (2) showed that 
(a = 0.8) and the exponent, b, was evaluated to be 0.4 from 
von Raman' s modification of Reynold's analogy for ML, = 1. The 
characteristic length in the Reynolds nuniber is the local diameter 
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which varies approximately with the boundary layer* One notable 
exception to this variation is in the entrance region of a nozzle 
where the boundary layer thickness may be small* 

The constant C was evaluated by determining the heat transfer 
coefficient at the throat of the nozzle shown in Figure 2* with 
the boundary layer analysis and solving for C in equation (l4) using 
this h . To insure that this equation would apply to other nozzle 

8 /VP - 1 


contours and conditions* a factor J , found from nozzle similarity 
studies in (2) was multiplied into the equation. If it is assumed 
that the specific heat* CL* and the Prandtl number* Pr* are constant 

Jr 

with temperature and that the gas properties p and ji are evaluated 
at a reference temperature to account for compressibility and/or heat 
transfer effects* equation (l4) can be expressed as 



This equation can be used to obtain a rapid estimation of the 
heat transfer distribution. It has been found to give good results 
except in the entrance region of nozzles with thin boundary layer 


thicknesses . 
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DESCRIPTION OF EXPERIMENTAL TECHNIQUE 

D iscussion of Experimental Apparatus 

The rocket engine used to obtain the experimental heat transfer 
data is shown in the schematic in Figure 3- The solid propellant was 
a composite of 83.3 percent by weight ammonium perchlorate ( oxidizer) , 
l4.4 percent polybutadiene acrylic acid (fuel) and 2.3 percent 
stabilizer. This propellant was selected because there are only small 
amounts of solid particles (carbon residue) in the exhaust products, 
which minimizes heat transfer by radiation, and because its combustion 
chamber temperature is relatively low (about 4500° F) . There is very 
little dissociation of the exhaust gases In this temperature range. 

Average chamber pressures of 220, 4io, and 742 psia were obtained 
by varying either the throat area or the volume of the combustion 
chamber and the amount of propellant. Table 1 gives the pertinent 
parameters for each of the chamber pressures used in the experiments. 

It was desirable to have a constant chamber pressure over the data 
taking period. To approximate this condition, the propellant was 
cast in a hollow cylinder configuration with the inner cylindrical 
surface being used as the exposed burning surface. Figures 4, 3 , and 
6 show the pressure traces for the three test conditions. It can be 
seen that a good approximation of constant chamber pressure with time 
was achieved. 

The convergent- divergent nozzle was made with a ZTA graphite 
insert and a steel housing. The internal surface contour was made smooth 
and continuous which allowed an assumption to be made that the 
beginning of the boundary layer coincided with the beginning of 
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convergent section of the nozzle. The ZTA graphite was chosen as the 
material for the insert because of its maehinability, resistance to 
erosion, and its thermal properties. Also by using the graphite 
insert, the experimental nozzle configuration was similar to nozzles 
used in full scale motors. 

Heat transfer measurements were made in the convergent region, at 
the throat, and in the divergent region of the nozzle for the three 
test conditions. A degree of redundancy was obtained by using two 
calorimeters diametrically opposed In the divergent region. These 
calorimeters also were used to determine whether the gas flow through 
the nozzle was concentric. The exact locations of all four calori- 
meters in the nozzle for each of the tests are shown in Figure 7. 

The calorimeters were made with a ZTA graphite core and a silica 
phenolic insulating sleeve. As can be seen in Figure 3, the graphite 
core was exposed directly to the flow of gases* By making both the 
insert and the calorimeter core out of the same material, there was 
very little disruption of the temperature distribution along the nozzle 
wall. The insulating sleeve was used to direct the heat flow along the 
longitudinal axis of the calorimeter. This allows the assumption of 
one- dimensional heat transfer to be made in the data reduction* The 
graphite core of the calorimeters was instrumented with five thermo- 
couples. The thermocouple nearest the exposed surface was composed 
of tungsten - 5 percent rhenium and tungsten - 2 6 percent rhenium wires 
5 mils in diameter. The other four thermocouples were made of 
platinum and platinum -13 percent rhodium wires also 5 oils in diameter. 
Figure 8 shows the locations of the thermocouples along the 
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longitudinal axis of the core. All calorimeters used in the tests 
were similarly constructed. 

By using five thermocouples in each calorimeter, a temperature 
distribution through the calorimeter could be obtained even if there 
were random failures of individual thermocouples. Also the thermal 
diffusivity of the graphite could be checked by specifying the tempera- 
ture-time history of one thermocouple in the data reduction analysis 
and comparing the calculated temperature distribution through the 
calorimeter with the distribution obtained from the experiment. This 
procedure will be demonstrated in a later section. 

Ideally, the surface temperature of the calorimeter should be 
directly measured. However* due to the extreme thermal environment 
at the surface, this measurement is very difficult to make. The 
surface temperatures of the calorimeters were determined by locating 
a thermocouple as closely as possible to the surface and extrapolating 

i 

the temperature data to the calorimeter exposed surface. 

All of the heat transfer measurements were taken from static 
firings of the rocket engine. Figures 9 and 10 show the engine in 
position on the thrust stand. The thermocouple data as well as chamber 
pressure aid thrust data were recorded on magnetic tape with a 
computer! 2 . ed data acquisition system. 

Thermal and Physical Properties 

The exhaust gas constituents were determined by assuming the 
products of combustion to be in chemical equilibrium at the prescribed 
pressure and enthalpy in the combustion chamber. Table 2 presents the 
conditions in the combustion chamber for the three tests and the 
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resulting mole fractions of each exhaust gas constituent as well as 
the specific heat ratio, y, and the specific heat at constant pres- 
sure, Cp, for the gas mirture. These values were determined with the 
use of a computer program described by Zeleznik and Gordon (21) The 
gas mixture was assumed not to vary in the nozzle (frozen flow) and 
local conditions at the test measurement locations were determined 
by expanding the mixture isentropically through the nozzle. 

The viscosity and thermal conductivity of the gas mixture were 

determined by using the equations suggested by Brokaw (6): 

v 



( 14 ) 



( 15 ) 



(16) 


(IT) 


These equations express the mixture viscosity and thermal conduc- 
tivity in terms of the viscosity and thermal conductivity of the 
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mixture costituents. The constituents that are starred in Table 2 
were accounted for in these calculations since they represent the 
major quantities in the mixture. The viscosity and thermal conductivity 
of each component were taken from data presented by Svebla (l6) . Their 
values were evaluated at the combustion chamber temperature. 

The properties of the ZTA graphite used in the calorimeter core 
that are required are the density, constant pressure, specific heat, 
and the thermal conductivity. The density and specific heat data 
were taken from ( IT ) - The thermal conductivity data were taken from 
the results of Wagner and Dauelsberg (l8) . Figures 11 and 12 show 
the variation of the specific heat and thermal conductivity of ZTA 
graphite with temperature. The thermal conductivity shown is for 
the "across the grain" direction of the graphite which was the direc- 
tion of heat flow in the calorimeters. It can be seen in the figure 
that the experimental data used to determine the curve is scattered. 

The thermal conductivity of any graphite is difficult to control. 

The effects of this variation on the test data will be discussed in 
the "Accuracy of Results" section. The density and specific heat of 
graphites do not vary nearly as much and can be controlled readily. 

Data Reduction 

All of the experimental data were recorded on magnetic tape so 
that the data could V ? reduced by computers , The thermocouple millivolt 
readings were converted into degree Fahrenheit according to National 
Bureau of Standards conversion tables. The temperature data were** * 
printed out in specified time increments over the firing time of the 
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rocket motors. The chamber pressure data were reduced in a similar 
manner . 

The heat flux into the calorimeter was determined by dividing 
the graphite core into 20 finite elements (Fig. 13), and specifying 
the temperature- time history of the third element. The temperature- 
time history was that of thermocouple 1, the thermocouple nearest the 
heated surface in each calorimeter. By using a finite element tech- 
nique , heat balance equations can be written for each of the elements 
over small time increments. For example , considering the heat flow 
to be one dimensional, the heat balance on element 1 is 




+ ft 


stored 


(16) 





+ pV, C 
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where is the temperature of element 1 calculated in the previous 

time interval. After the heat balance equations are written for each 
of the elements, these equations may be solved simultaneously for the 
heat transfer coefficient. 

A computer program was used to solve for the heat transfer 
coefficient versus firing time for each of the calorimeters. The 
program utilizes the Gauss-Jordan method to solve the simultaneous 
linear heat balance equations. In addition to a temperature- time 
history of one of the elements, necessary input data include the 
physical dimensions of each element and the thermal properties of the 
material (K, p, and C ) which may vary with temperature. Output 

Jr 

includes the heat transfer coefficient, the heating rate at che surface 
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of the calorimeter, and the temperature gradient along the longi- 
tudinal axis of the calorimeter for each specified time increment. 

Radiative heating from the exhaust gases was not accounted for in 
the analysis. This was based on the fact, as discussed by Kuby (ll), 
that gases radiate energy in finite frequency bands as opposed to most 
solid bodies which radiate in a energy continuum and therefore the 
total integrated value of emitted energy is much less for gaseous 
radiation. The propellant used in the tests was chosen because the 
amount of solid particles in its exhaust gases is small enough to be 
neglected. By not accounting for the radiation from the exhaust 
products, the data reduction is simplified, and the results are not 


effected substantially 
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RESULTS 


Presentation and Discussion of Results 
Experimental heat transfer coefficients determined for average 
chamber pressures of 220, blO, and Jk2 psia are presented. The experi- 
mental coefficients are compared with theoretical coefficients deter- 
mined from the two methods of Bartz. In determining the theoretical 
coefficients with the simultaneous solution of the boundary layer 
momentum and energy equations, the skin friction coefficient was 
evaluated at the arithmetic mean between the free stream and wall 
temperature as well as at the free stream temperature. Also, the 
theoretical coefficient was evaluated with an "interaction component" 
value of 0. 1 as recommended by Bartz and also with a value of zero 
which essenti all y decouples the momentum and energy equations. 

Other data presented are the calculated temperature distribution 
derived from the temperature data from the thermocouple 1 of each 
calorimeter. This calculated temperature distribution is compared with 
the temperature readings of the thermocouple mounted along the longi- 
tudinal axis of the calorimeters to determine the accuracy of the 
experimental data. The deviations between thermocouple readings and 
the calculated temperatures are discussed in the "Accuracy of Results" 
section. The experimental heating rates are also presented. 

The thermocouples are numbered consecutively from one to five with 
thermocouple 1 located nearest the heated surface of each calorimeter. 


Average Chamber Pressure of 220 psia 

Experimental heat transfer data were obtained at a local area 

% 

to throat area ratio , •— , cf 1.785 i n the convergent section of the 
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A 1 

nozzle , at the nozzle throat, -r— ® 1, and In the divergent section at 
A 1 ^ 

— = 2 . 569 , Figures ik through 22 show the temperature response of the 
A t 

three thermocouples nearest the heated surface In each of the calori- 
meters, the experimental heating rates, and the experimental heat 
transfer coefficients. The oscillatory nature of the heating rate and 
heat transfer coefficient curves is due to the data reduction technique 
rather than physical conditions in the test. 

The heating rate curves show a decline in the magnitude of the 
heat flux with increasing time. Tills is to be expected since the 


driving potential 


’ ( T a» - \) ■ 


becomes smaller with increasing time. 


However, the heat transfer coefficients calculated from the measured 
temperatures in the convergent section and at the throat also decreased 
with time. In Figure 25 it can be seen that both of Bartz’s methods 
overpredict the heat transfer coefficient throughout the nozzle. The 
heat transfer data from the latter part of the firing can be correlated 
by using a laminar flow equation used in (12) : 



I *f0*5 
3 *e D 


K 


.-0.067 

Fr 


(18) 


This is Pohlhausen’s equation for laminar flow where the characteristic 
length has been changed to the local diameter since the flow Is 
internal. The heat transfer coefficient determined from the experi- 
mental data in the divergent region of the nozzle did not vary as much 
as the throat and convergent heat transfer coefficients. The predic- 
tions from the simultaneous solutions, where the shin friction 
coefficient was evaluated at the free stream temperature, and from the 
laminar flow equation bracketed the test data in the divergent region. 
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The leveling of the slope of the temperature curves of thermo- 
couple 1 in the convergent and throat calorimeters seemed unusual. In 
order to determine if the temperature data were valid, a temperature- 
time response of thermocouple 1 in each calorimeter was assumed as shown 
in Figures 24 and 25 . These assumed temperature responses were used 
to compute the temperature response throughout the calorimeter with the 
use of the finite element computer program. Figures 24 and 25 show that 
on comparing the results using the assumed response with thermocouples 
2 and 5 in the convergent calorimeter and thermocouple 2 in the throat 
calorimeter, the assumed temperature response was incorrect. The 
calculated responses were higher than the experimental responses at 
the locations where temperature were measured. These comparisons 
lend credance to the experimental data which show that the heat transfer 
coefficient becomes smaller with time. 

It is known (4) that cooling the laminar boundary layer increases 
e stability, i.e., increases the critical Reynolds number for 
transition to turbulent flow. The boundary layer in this test was 
cooled by the nozzle wall throughout the test but to a lesser extent 
as time increased. Thus, it could be assumed that the conditions for 
laminar flow became less suitable at later times in the firing. The 
experimental data, however, indicate a tendency from turbulent or 
transition flow to laminar flow with increasing time. 

It is generally assumed (4) that boundary layer heat transfer coef- 
ficients are affected by wall temperature. This is accounted for by 
evaluating the skin friction coefficient and the gas transport properties 
at some intermediate temperature between the adiabatic and wall tempera- 
ture. Figure 23 shows the difference between evaluating the skin friction 
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coefficient at the free stream temperature and the arithmetic mean 
of the free stream and wall temperature for turbulent flow. Increasing 
wall temperature would tend to lower the heat transfer coefficient but 
not to the extent experienced in the test, 

A favorable pressure gradient tends to stabilize a laminar boundary- 
layer and since the local pressure decreases along the length of the 
nozzle, a turbulent boundary layer could revert to laminar flow. In 
Figure 2 6 the experimental data from the 220 psia chamber pressure 
test is compared with experimental data presented in (12). The figure 
is a plot of the parameter, Stanton number multiplied by the Prandtl 
number to the o£ power versus free stream Reynolds number based on the 
local diameter. It can be seen that the data reported on herein falls 
within the transition region as defined in (12). These results suggest 
that it is possible for laminar or transition flow to occur in a rocket 
motor nozzle even though the exhaust gases originate from a combustive 
process within the motor chamber that is highly turbulent. 

Average Chamber Pressure of 4l0 psia 

The nozzle used in this test was the same as that used in the 
220 psia test with measurements made at the same locations. The addi- 
tional pressure was generated by coupling two motor cases together 
with the head end case containing only half as much propellant as the 
second case. This method gave a chamber pressure that varied more than 
in the 220 psia test. 

Figures 27 through 35 show the temperature response of the three 
thermocouples nearest the surface in each of the calorimeters, the 
heating rates, and the experimental heat transfer coefficients. 
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Figure 36 shows a comparison of the experimentally determined coef- 
ficients with the several methods of predicting the coefficients. It 
can be seen that, as in the 220 psia test, both of the Bartz methods 
overpredict the heat transfer coefficients, with the technique of 
evaluating the skin friction coefficient at the free stream temperature 
coming nearest to correlating the data. The heat transfer coefficient 
at the throat of the nozzle decreased sharply with time. This is also 
occurred to a lesser extent at the convergent calorimeter. The local 

Reynolds number at the convergent, throat, and divergent calorimeter 
6 6 6 

was O.T23 X 10 , 1.02 x 10°, and 0 .793 x 10, respectively, indicating 
turbulent flow when compared to the data in (12). The drastic reduc- 
tion of the heat transfer coefficient at the throat is questionable 
in view of the fact that the Reynolds number at the throat indicates 
turbulent flow and that the magnitude of the coefficient drops slightly 
below that of the convergent calorimeter. The initial 1.5 to 2.0 seconds 
of data do indicate that the Bartz methods slightly overpredicted the 
results of the test at the throat. 

There was very little variation with tine in the heat transfer 
coefficient on the divergent calorimeter. By evaluating the skin 
friction coefficient at the free stream temperature, the heat transfer 
coefficient was very nearly predicted. The convergent calorimeter 
data did show some decrease in the heat transfer coefficient. The 
heat transfer data at this point was slightly overpredicted by using 
the free stream skin friction coefficient in the simultaneous solution. 

As mentioned before, there was some variation in the chamber 
pressure. An average chamber pressure of ^10 psia was used in the 
analytical solutions. The chamber pressure varied from Wo psia in the 
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initial stages of the firing to $60 psia just before burnout. Using 
the Nusselt number correlation equation, where the heat transfer 
coefficient is proportional to the 0.8 power of the chamber pressure, 
to determine the effect on the heat transfer coefficient, it is found 
that the chamber pressure variation would theoretically cause the 
heat transfer coefficient to vary 6.5 percent above to 9-8 percent 
below that for an average chamber pressure of 4l0 psia. The change 
in chamber pressure could partially explain the reduction in heat 
transfer coefficients at the convergent calorimeter and at the throat 
calorimeter. 


Average Chamber Pressure of 7k2 psia 
The nozzle contour for this test was slightly different for this 
test in order to increase the chamber pressure to 7^2 psia. The 
contour still included at 45° half angle convergent region and a 15° 
half angle exit cone, but had a smaller throat diameter. Heat transfer 
measurements were made at an area ratio, — = 2.067 in the convergent 
region, -r— = 1, at the throat, and at — ^ 3*764, in the divergent 
region. Figures 37 through 45 show the temperature responses of the 
three thermocouples nearest the heated surfaces in each of the calori- 
meters, the heat fluxes calculated from the temperature responses, and 
the resulting heat transfer coefficients. Figure 46 compares the experi- 
mental results with the various analytical prediction methods. This 
figure shows the data to be slightly overpredicted when the free stream 
skin friction coefficient is used in the simultaneous solution. The 
convergent region dat * showed very little variation of the heat transfer 
coefficient during the burning of the motor although the overprediction 
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was the largest at this measurement location. The heat transfer coef- 
ficient at the throat was in the range of the predictions that used the 
free stream skin friction coefficient initially and decreased as the 
burning progressed. The divergent calorimeter data agreed very nearly 
with the predictions of the method stated above. 

Accuracy of Results 

The accuracy of the experimental results can be broken down into 
two parts: 

(1) Accuracy of the measured data. 

(2) Accuracy of the method of data reduction. 

Factors which affect the accuracy of the measured data are the ability 
to determine the true location of the thermocouples with respect to the 
nozzle internal wall and the error introduced by the recording equipment. 
Factors which affect the accuracy of the data reduction are how well 
the properties of the ZTA graphite can be determined, i,e., the density 
and the constant pressure specific heat and thermal conductivity as a 
function temperature, and how accurate the method of data reduction is. 

The location of the thermocouples in the calorimeter could be 
determined to x^tthin ±.003 inch by using close tolerances in the fabri- 
cation of the calorimeter and X-raying each calorimeter. Figure 47 is 
a typical X-r*y photograph. The tolerances on the installation of 
the calorimeter in the nozzle were such that the true location of the 
thermocouples with respect to the nozzle internal wall could be deter- 
mined to within ±.008 inch. This resulted in a band of uncertainty of 
the temperature distribution along the length of the calorimeter as 
shown in Figure 48 . 
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The thermocouples used in the tests could measure true temperature 
within 1 percent. The noise in the recording channels of the measuring 
equipment could be filtered so that there was negligible error in the 
recording of the temperatures. 

As stated before, the experimental data on the thermal conductivity 
of ZTA graphite is scattered (Fig. 12). In attempting to fit the 
temperature data from the thermocouple 1 with the experimental 
temperature distribution, it was necessary to vary the thermal conduc- 
tivity of the graphite. The curves shown in the figures that present 
the temperature data versus time represent a "best fit." Since the 
thermal conductivity of the graphite varies with temperature, the 
thermal conductivity versus temperature curve had to be varied rather 
than a thermal conductivity value. The density and constant pressure 
specific heat versus temperature of ZTA graphite are accurately 
know (l8) and error in their values were not considered in the analysis. 

In order to determine the effect of the ±.008 inch error in the 
location of the thermocouple 1, an error analysis was made by 
varying its location in the data reduction* The calculated heat input 
to the nozzle wall was found to vary ±5 percent during the early times 
in the test with the variation decreasing to ±2 percent during the 
later times in the test. This variation was in the calculated heat 
input only. As previously stated, the experimental data was not 
perfectly fitted with the calculated data based on the temperature data 
of thermocouple 1, Figures 48 through 50 show the percent difference 
between the heat input based on the temperature data of thermocouples 
1, 2, and 5 and the calculated heat input for late times in the tests. 
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The percent difference was determined by comparing the slopes of 
curves at the heated surface of the calorimeter since the heat input 


is: 


% =* -K 


ST 
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surface 


( 19 ) 


The thermal conductivity was not varied since the difference in surface 
temperature of the two curves was not large. 

In the figures that show the calculated heat input and the 
calculated heat transfer coefficients , there is oscillation in the 
curves caused by the data reducing process rather than any physical 
phenomena occurring in the tests. This was caused by the large change 
in temperature of the graphite in a short time, the temperature varia- 
tion of the specific heat and thermal conductivity of the graphite, and 
the fact that a finite difference technique was used. 

In summary, the heat transfer results from the tests were not 
exact results. However, the data did show trends which were presented 
and discussed previously. 
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SUMMARY AND CONCLUSIONS 

The experimental heat transfer coefficients obtained by firing a 
solid propellant rocket motor at chamber pressures of 220, 4l0, and 
7^-2 psia were generally below the predictions made by the Bartz methods 
of solving the boundary layer energy and momentum equations simultaneously 
and the Russelt number correlation equation. 

The test at 220 psia chamber pressure demonstrated that it is 
possible to obtain laminar or transitional boundary layer flow in a 
solid propellant motor even though a turbulent combustive process is 
occurring in the motor chamber. The best agreement was found in the 
divergent region of the nozzles at all three pressures where the experi- 
mental coefficients were only slightly below the predictions from the 
simultaneous solution and evaluating the skin friction coefficient at 
the local free stream temperature. 

Data from the convergent and throat calorimeters indicated the 
heat transfer coefficient decreases with increasing wall temperature. 

This variation was more pronounced than the variation in the predic- 
tions when the skin friction coefficient was evaluated at the arithmetic 
mean of the free stream and wall temperature (film temperature). The 
predictions using the film temperature skin friction coefficient were 
always above the experimental heat transfer coefficients in magnitude. 
Predictions using the skin friction coefficient evaluated at the free 
stream temperature were in the same range of the experimental data, 
but they do not account for wall temperature variation. 

These results tend to substantiate the data in the literature in 
the respect that predicting heat transfer rates in the inlet portion of 
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a rocket motor nozzle is difficult (k, 19, 20 ) and that transitional 
and laminar flow along the nozzle wall is possible (l, 5, 12) . The 
data in (5* 12) were collected from nozzles similar to the nozzles 
used in the tests reported on here and for similar combustion chamber 
conditions (F^ on the order of 200 psia) . The data from (l) showed 
the transitional flow phenomena to occur at a lower combustion chamber 
pressure (75 psia) . 

These test data indicate that by using either the Bartz simultaneous 
solution or the Uusselt number correlation equation recommended by 
Bartz, calculations of the heating load along a nozzle wall would be 
conservative in that the analytical heat transfer coefficient is higher 
than the experimental coefficient. It is felt, however, that the degree 
of conservatism will hamper efforts to have an efficient nozzle design. 
Further analytical and experimental efforts are needed to refine and 
improve Bartz 1 s methods. 

It should be noted that these sets of data are not a complete test 
of the Bartz methods since the data was derived from the use of one 
type of nozzle. A complete conclusion may be drawn when test data are 
available from a variety of propellants and nozzles )ised on solid 
propellant motors. 
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Table 1. Varlati^r of Motor Characteristics to Achieve Different Chamber Pressures 
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Table 2 


Exhaust Gas Properties 


Average chamber 
pressure, psia 


Chamber gas 
temperature, °R 


Specific heat 
ratio 


Specific heat 
constant. Pressure 
Btu 


Viscosity 


Thermal conductivity 
Btu ft 


Exhaust gas 
constituents, mole 
fractions 


Mole fractions of constituents did not vary with chamber pressure 


Note 





Dimensions in Inches 


glow Conditions 

p Q = 300 paia C p * 0.567 Btu/lb °E 7 u 1.2 

T 0 = 4,500° H N pr = O .83 H ~ T 0,65 

n 0 - 1.3 X 10" 8 lb sec /ft 2 

Figure 8. Nozzle for evaluation of constant in heat-transfer coeffi 
oient equation (from ref. 4) 



To pressure 
transducer 



Figure 3- Experimental solid -propellant motor schematic 



Time, see 


Figure K Chamber pressure vs time - P avg 






. 15 . 



220 psia 


Tim e, sec 


Figure % Chamber pressure vs time - F avg 



klO psia 









o rv 


Throat 
Diameter* in 


Calorimeter locations within nozzle 
















Kote; All thermocouples were located on graphite core centerline 


Dimensions for 
convergent and 
calorimeters 



Five thermocouples 
eq ua lly spaced 0#28 apart 


Five thermocouples 
equally spaced 0,253 apart 



Silica phenolic 
sleeve with 
0,10 wall 

ZTA graphite 
core 


/ 


All dimensions in inches 


Dimensions for 

divergent 

calorimeters 

> 


Figure 8. Calorimeter dimensions 



Figure 9 


Experimental solid -propellant motor mounted in thrust stand 










Convergent and throat 
calorimeters 


Element Thickness, in. 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


.0210 

.0210 

.0360 

.0475 

.0625 

.0625 

.0625 

.0625 

.0625 

.0625 

.0625 

.0310 

.0625 

.0940 

.0625 

.0625 

.0625 

.1040 

.1040 

.1040 


Divergent 

calorimeters 


Heat in 



Side and end 
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Figure 13. ZTA graphite core finite elements 
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Figure Ik. Comparison of experimental and calculated temperature response of thermocouples 1 
2, and 3* Convergent calorimeter - P c = 220 psia 
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Figure 17 . Comparison of experimental and calculated temperature response of the? locouples 1 
and 2. Throat calorimeter - P r - 220 psia 
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Figure 24. Extrapolated thermocouple number 1 response. Convergent calorimeter - P c - 220 psia 
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Figure 25 . Extrapolated thermocouple number 1 response. Throat calorimeter - ? c = 220 psia 






Figure 27- Comparison of experimental and calculated temperature response of thermocouples 1 
2, and 3* Convergent calorimeter - P c = 410 psia 
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Figure 29* Convective heat -transfer coefficient vs burn time. Convergent calorimeter 
P c = klO psia 
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Figure Jl. Heating rate vs burn time. Throat calorimeter - P c = 410 psia 
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Figure 33 * Comparison of experimental and calculated, temperature response of thermocouples 1 
2, and 3 * Divergent calorimeter - P c = klO psia 
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Figure k-0. Comparison of experimental and calculated temperature response of thermocouples 1 
2. and 3- Throat calorimeter - P = 7 42 psia 













Figure *4-5. Convective heat-transfer coefficient vs burn time. Divergent calorimeter 
P c = 7*4-2 psia 
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